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Abstract 



We study the phase diagram and critical behavior of a one dimensional 
three species monomer-monomer catalytic surface reaction model. Static 
Monte Carlo simulations are used to roughly map out the phase diagram 
consisting of a reactive steady state bordered by three equivalent unreac- 
tive phases where the surface is saturated with one monomer species. The 
transitions from the reactive phase are all continuous, while the transitions 
between poisoned phases are first-order. Of particular interest are the bi- 
critical points where the reactive phase simultaneously meets two poisoned 
phases. A mean-field cluster analysis fails to predict all of the qualitative 
features of the phase diagram unless correlations up to triplets of adjacent 
sites are included. Scaling properties of the continuous transitions and the 
bicritical points are studied using dynamic Monte Carlo simulations. The 
transition from the reactive to a saturated phase shows directed percolation 
critical behavior, while the universal behavior at the bicritical point is in the 
even branching annihilating random walk class. The crossover from bicritical 
to critical behavior is also studied. 
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I. INTRODUCTION 



Nonequilibrium models with many degrees of freedom whose dynamics violate detailed 
balance arise in studies of biological populations, chemical reactions such as heterogeneous 
catalysis, fluid turbulence, and elsewhere. The macroscopic behavior of these models can 
be much richer than that of systems in thermal equilibrium, showing organized macroscopic 
spatial and temporal structures like pulses or waves, and even spatiotemporal chaos. Even 
the steady state behavior can be far more complicated, involving for example scale invari- 
ance at generic parameter values, and critical behavior distinct from any equilibrium models. 
However, like their equilibrium cousins, systems at continuous transitions between nonequi- 
librium steady states show universal behavior that is insensitive to microscopic details and 
depends only on properties such as symmetries and conservation laws. 

One place where such nonequilibrium models appear is in the study of chemical reac- 
tions occurring on catalytic surfaces, which show a variety of interesting behavior including 
nonequilibrium phase transitions, temporal oscillations, spiral waves, and chemical chaos [|I| . 
In order to help understand these complicated processes, a number of simple models have 
recently been proposed that attempt to capture the essential physics 0. 

Ziff, Gulari, and Barshad (ZGB) proposed a monomer-dimer reaction model to explain 
some features of CO oxidation on a noble metal surface [Q. In their model, monomers 
representing CO molecules and dimers representing O2 molecules adsorb on a lattice. Im- 
mediately upon adsorption, the O2 dimers dissociate into two O monomers. CO monomers 
and O monomers occupying nearest neighbor sites then react to form a CO2 molecule that 
immediately desorbs, leaving two vacant lattice sites. In the limit of infinitely fast reactions 
(the adsorption controlled limit), where the only parameter of the model is the relative 
adsorption rate of CO molecules ycoi they found in two dimensions that there are three 
phases: An O2, or dimer poisoned state for yco < a CO, or monomer poisoned phase 
for yco > 1/2, and a reactive phase for yi < yco < At yi the fraction of each species 
changes continuously, indicating that the dimer poisoning transition is continuous. At 1/2 
the monomer poisoning transition is first order, with the densities of the different species 
changing discontinuously. In one dimension, the ZGB monomer-dimer reaction model has no 
reactive phase, only monomer-poisoned and dimer-poisoned phases separated by a first-order 
transition [Q. 

An even simpler catalytic reaction model can be constructed by replacing the dimer 
species in the ZGB model with a second monomer species. This monomer-monomer model 
has a long history , and in fact certain analytic results for this model have been obtained in 
the reaction-controlled limit of the model [Q. In this model two different monomer species, 
call them A and B, adsorb on a lattice where nearest neighbor AB pairs react and an AB 
molecule desorbs. However, the phase diagram for this model does not contain a reactive 
steady state in any number of dimensions, neither in the adsorption controlled nor in the 
reaction controlled limit. The phase diagram consists only of A and B poisoned states, and 
a first-order transition between them. 

The dimer poisoning transition in the ZGB model is one of the most common types 
of continuous phase transitions in nonequilibrium models. It is a transition to a single 
absorbing, noiseless, steady state, the term absorbing indicating the state cannot be left 
once it is reached. Other examples include directed percolation (DP) the contact 
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process 0, auto-catalytic reaction models and branching annihilating random walks 



with odd numbers of offspring [|TT],|12|- Both renormalization group calculations |^,13| and 
Monte Carlo simulations P-p!^,p!^ show that these models form a single universality class 
for a purely nonequilibrium model with no internal symmetry in the order parameter. 

Recently, a number of models with continuous adsorbing transitions in a universality class 
distinct from directed percolation have been studied. These models include probabilistic 
cellular automata models studied by Grassberger et al. []T3[, certain kinetic Ising models 
||16| , the interacting monomer-dimer model ? and branching annihilating random walks 

with an even number of offspring (BAWe) [TT.TI]]. All of these models except for the BAWe 
have two equivalent absorbing states indicating the importance of symmetry of the adsorbing 
state to the universality class. However, the universal behavior of this new class is apparently 
controlled by a dynamical conservation law. If the important dynamical variables in this class 
are defects represented by the walkers in the BAWe model and the walls between different 
saturated domains in the other models, the models have a "defect parity" conservation law 



15 where the number of defects is conserved modulo 2. Recent field theoretic work confirms 



this viewpoint [pO| ]. 
In a recent Letter 
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we introduced a monomer-monomer reaction model with three 
This model could represent either a system with three different 



different monomer species, 
chemical species or an auto-catalytic reaction system in which one chemical species can 
adsorb on three different types of surface sites. Using static and dynamic Monte Carlo 
simulations, we determined the phase diagram and studied the phase transitions in the one 
dimensional version of the model, and showed that it has continuous adsorbing transitions 
to both one and two equivalent noiseless states. It is therefore is a good model to study the 
role of symmetry in adsorbing phase transitions. 

In this paper we expand those results, providing more details of our simulation methods 
and of the results, again restricting our consideration the the one-dimensional version of the 
model. We also include a mean-field cluster analysis of the model including up to triplets of 
adjacent sites. The paper is organized as follows. In the next section we define the model and 
show the phase diagram of the model, as determined by simulations. The following section 
presents the mean-field analysis. Section IV contains the details and results a detailed 
Monte Carlo study of the dynamic scaling behavior at the various phase transitions, and of 
the crossover behavior between the different types of scaling behavior. In the last section 
we summarize our results. 



II. THE MODEL 

Our three species monomer-monomer model is defined by two fundamental dynamic 
processes: (a) monomer adsorption at sites of a substrate, and (b) the annihilation reaction 
of two dissimilar monomers adsorbed on nearest-neighbor sites of the substrate. Here we 
consider the model only in the adsorption controlled limit where process (b) occurs instan- 
taneously. Calling the monomer species A, B and C, the parameters in the model are then 
the relative adsorption rates of the different monomer species pa, Pb, and pc, such that 
Pa + Pb + Pc = 1- Using static Monte Carlo simulations to get the rough picture, and 
refining it with dynamical Monte Carlo studies described below, we find the ternary phase 
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diagram for the model is shown in Fig. |l[ In this figure, the horizontal axis corresponds 
to the relative adsorption rate of A and B monomers pab = Pa/{pa +Pb)- The absorbing 
phases, where one monomer species saturates the chain, occupy the corners of the phase 
diagram. In the center of the phase diagram is a a reactive steady state. There are contin- 
uous phase transitions from the reactive phase to the saturated phases, but the monomer 
densities undergo discontinuous, first-order, transitions from one saturated state to another. 
The points where the reactive phase and two saturated phases meet are bicritical points 
where two lines of continuous transitions meet a line of first-order transitions. 



III. MEAN-FIELD THEORY 

To analyze the kinetics of the three species monomer-monomer model, it is useful to 
perform a mean-field analysis. While such analysis neglects long-range correlations and thus 
cannot be expected to properly predict critical properties, it should properly predict the 
qualitative structure of the phase diagram, including the existence of continuous transitions 
and multi-critical points. The mean-field analysis also provides a starting point for studying 
the importance of such fluctuations, which, of course, become particularly important near 
continuous phase transitions. The mean-field approach we use studies the time evolution 



of clusters of sites, the approximation coming in truncating the probabilities of observing 
clusters of larger size into probabilities for smaller size clusters. The simplest form is the site 
approximation where probabilities of observing certain nearest neighbor pairs is replaced by 
the produce of the average site densities. Better approximations can be obtained systemat- 
ically by replacing the actual configuration of larger clusters, i.e. pairs, then triplets, and 
so on, with the average density of those clusters. The analysis presented below includes 
clusters consisting of up to triplets of adjacent sites. 



A. Site approximation 

At a particular time, a lattice with N sites will have Ny vacancies, the remaining sites 
being filled with N^, Nb, and A^^^ numbers of A, B, and C monomers respectively. The 
density of A monomers is xa = (Na/N), with corresponding definitions for 5, C and vacant 
{V) sites. We have the obvious constraint 

Xv + Xa + xb + xc = I- (1) 

In the site approximation all correlations are neglected, so that Xy is the probability that 
a given pair of lattice sites are occupied by two vacancies. The rate equations for the A 
monomer density is 



dx 



^ - PaXv{^ - Xb- xcf - {pb + Pc)xv 1 - (1 - xaY (2) 



dt 



with similar equations for Xb and Xq- The first term on the right hand side of Eq. is the 
rate of A monomer adsorption multiplied by the probability that an adsorbing A monomer 
will find a vacant site that has no S or C monomers adsorbed on adjacent sites. The second 
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term is the rate at which B or C monomers find a vacant site with at least one adjacent 
adsorbed A monomer to react with. 

Equations (0) have steady state solutions corresponding to each of the three poisoned 
states, as well as one corresponding to the reactive steady state. To find the site approxi- 
mation phase diagram, we analyzed the stability of those solutions as a function of the rates 
{Pa} by examining the eigenvalues of the Jacobian matrix for linearized rate equations. 

For example, the Jacobian matrix for the A poisoned state has two zero eigenvalues and 
one eigenvalue of Pb+Pc~Pa = (1 ~ 2p^). This third eigenvalue shows that the A poisoned 
state is stable only for pa > 1/2. Corresponding results hold for the other poisoned states, 
leading to the site approximation phase diagram shown in Fig. 0. 

As the phase boundaries are approached from the reactive phase, the monomer densities 
vanish continuously, indicating a continuous transition to an absorbing state. The points on 
the edge of the phase diagram where two different poisoned phases meet the reactive phase 
are bicritical points. 



B. Pair approximation 

We improve the site approximation by properly accounting for the correlation of nearest 
neighbor pairs and approximating the correlations of triples and larger clusters. We define 
Nij as the number of bonds connecting nearest neighbor sites occupied by i and j monomers 
{A, B, C, or V), where the monomer i occupies the site to the left of monomer j, and 
we have Nij = Nji. Since we are studying one dimension, the number of bonds equals the 
number of sites A^, so the bond densities are defined by 

AT.. + A^.. 



and 



Xii- ^ 



There are seven different allowed types of bonds: V-V, A-A, B-B, C-C, A-V, B-V, and C-V. 
Other types of bonds, A-B, A-C, and B-C, are forbidden in the adsorption controlled limit 
we are considering. The densities satisfy the constraint 

^VV + ^AA + ^BB + ^CC + ^AY + ^BV + ^CV = 1, (3) 

SO only six of the are independent. The A monomer density is given by xa = xaa + \xav-i 
with similar expressions for the B and C densities. 

To determine the equations of motion of the pair densities it is useful to distinguish 
between the different types of events that change the configuration. For example, if an A 
monomer attempts to occupy a site, it can (1) stick, (2) react with a i? or (3) react with 
a C, which we indicate respectively with the shorthand : (1) A |, (2) A [ AB |, and (3) 
A [ AC t- The rate equations can be written as 



dt ^ 
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where a refers to the event type, and Ax •J'' is the change in ij bond density arising from an 
event of type a. 

To find the different bond density changes note that the probabihty P{i\j) for a site to 
be occupied by a monomer (or vacancy) of type i, given that one of its nearest neighbors is 
of type j, is 

pm = w = ?- 

Mi 2xj 

for i j, and 

pm = ^ 



The various Axj^^ are given in Table 1, where 



z,v = Pii\V) + PiV\V) 



Xvv + l^iV 



Xv 



is the probabihty that the site to the left of a vacant site is occupied by either an i type 
monomer or a V. The density changes due to the other event types are found by permutation. 
Thus, the rate equations are 

dXAA XaaXav r /I , n , /i , m 

__ = xav zav 2x^ ^^^^ + + zcv)\ 

dx x^ 

= Pa zavC^xvv - Xav) - [Pb{^ + zbv) +Pc(1 + zcv)] (4) 

The other equations can be found by permutation, except for xyy which can be found 
using Eq. (0). 

Multiple steady state solutions to the set of six coupled bond density rate Eqs. (0) cor- 
respond to the reactive state (which can be found numerically), as well the three poisoned 
states. In principle, to find the phase diagram a stability analysis of those steady state solu- 
tions could be performed. However, we instead simply solved the six equations numerically 
as a function of the parameters pab and pc, and looked for the transitions to the poisoned 
states. The results are shown in Fig. |[ The densities of the different monomer species still 
change continuously as the phase boundaries are approached, indicating that the transitions 
are continuous. While the phase boundaries are now curved as they are in the actual phase 
diagram, the bicritical points are still on the edge of the phase diagram unlike the actual 
phase diagram. 



C. Triplet approximation 

The mean-field theory can be refined even further by considering larger clusters. However, 
this systematic process rapidly increases in difficulty. But since even the pair approximation 
failed to predict that the bicritical points occur on the interior of the phase diagram, we 
pushed the cluster expansion one step further and analyzed the model in the triplet ap- 
proximation. In this approximation, clusters of three adjacent sites are considered, thereby 
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including the effects of correlations up to that level. The details of the calculation are pre- 
sented in the Appendix, but here we summarize the results. In one dimension, there are 19 
different allowed triplets. However, 4 different constraints reduce the number of independent 
triplets to 15. Numerically solving the rate equations for the densities of those 15 different 
triplets simultaneously, we find solutions corresponding to the reactive steady state, as well 
as the poisoned states. The phase diagram, calculated as for pair approximation, is shown 
in Fig. ^. Finally at this level of approximation all of the qualitative features of the actual 
phase diagram are predicted. In particular, the bicritical points appear on the interior of the 
phase diagram and there are first order lines between the poisoned phases. However, note 
that the size of the poisoned phases is still underestimated by the mean-field cluster analysis, 
even in the triplet approximation. For example, the bicritical point on the pab = 0.5 line 
occurs at about pc = 0.02 in the triplet approximation, whereas in actuality it occurs at 
about Pc = 0.12. This indicates that fluctuations, which are still not fully accounted for in 
mean-field theory, stabilize the poisoned phases. 

IV. SIMULATIONS 

To further investigate the three species monomer-monomer model we also used time- 
dependent Monte Carlo simulations. Particularly useful for studying critical properties, the 
method is a form of "epidemic" analysis |jl2|,|19|,^ in which the average time evolution of 
a particular configuration that is very close to an adsorbing state (defect dynamics), or 
very close to a minimal width interface between two different adsorbing states (interface 
dynamics), is measured by simulating a large number of independent realizations. Using 
this technique we determined the universality classes of the critical and bicritical points, and 
studied the critical dynamics of interfaces between the two symmetric saturated states at the 
bicritical points, and the crossover from bicritical to critical behavior, including measuring 
the crossover exponent 0, as well as the subcritical behavior at the first-order lines. 

Because monomers can adsorb only at vacant sites, and the total number of vacancies 
on the lattice is usually very small, instead of randomly picking a site to attempt to adsorb 
on, it is much more efficient to use a variable time algorithm in which the adsorption site is 
randomly picked from a list of vacant sites. The species of monomer chosen for adsorption is 
then randomly picked according to the relative adsorption rates {pa}, and the time length of 
a step is l/nv{t) where nv{t) is the total number of vacancies at that time. Thus, on average 
there is one attempted adsorption per lattice site per unit time. We always start with a 
lattice big enough that the active region will never reach the boundaries; it is effectively an 
infinite lattice. 

During the simulations we measured the survival probability P{t), defined as the proba- 
bility that the system had not poisoned by time t, the average number of vacancies per run 
{nv{t)), and the average mean-square size of the active region per surviving run {R'^{t)). At 
a continuous phase transition as t — > cxd these dynamical quantities obey power law behavior 

P(t)~r^ (r^v(^))~^^ {R\t))r^t\ (5) 

Plots of the logarithms of these quantities versus the logarithm of the time, such those shown 
in Fig. 1^ yield a straight line at the phase transition, and show curvature away from the 
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transition. 

Precise estimates of the location of the critical point and of the exponents can be made 
by examining the local slopes of the curves on a log-log plot. The effective exponent 6{t) is 
defined as 



= {In [P(t)/P(t/6)]/ In 6}, (6) 

with similar expressions for ri{t) and z{t). At the critical point, a graph of the local slope 
versus should extrapolate to the critical exponent, with a correction that is expected ||2^ 
to be linear in t~^. Away from the critical point, the local slope curve should show strong 
curvature away from the critical point value as — > 0. 



A. Critical dynamics 



Figure |^ shows the data for the three dynamic quantities near the phase transition to 
the C saturated phase at pab = 0.5 plotted against time on a log-log scale. This data was 
calculated from 10^ independent runs of up to 10^ time steps at each parameter value. As 
expected, right at the critical point the line is straight, indicating power law scaling, and 
away from the critical point the lines show curvature. 

The exponents and the location of the critical point are easily and precisely determined by 
taking the local slopes of this data, which are shown in Fig. ^. We find a critical C monomer 
adsorption rate of pc = 0.39575(10), and that the critical exponents are 6 = 0.16(1), r] = 
0.31(1), and z = 1.255(15). These values are consistent with our expectation that this 
transition should be in the DP universality class, for which the exponents are 6 = 0.1596(4), 
rj = 0.3137(10), and z = 1.2660(14) [^]. We found similar exponents for the adsorbing 
transition at a number of other points along the lines separating the reactive phase and 
the saturated states (see the discussion of the crossover from bicritical to critical behavior 
below), indicating the transition between the reactive phase and any single saturated phase 
is always in the DP universality class. 



B. Bicritical defect dynamics 



The same kind of analysis at the bicritical point at p^B = 0.5, using an initial condition 
of a vacancy in an A-saturated phase, yields a bicritical point at Pc = Pc — 0.122(1). The 
exponents at the bicritical point are very different, which we expect given the presence of 
two-symmetry equivalent saturated phases. From 5 x 10^ runs of up to 10^ time steps we 
found the local slope data shown in Fig. |^, yielding values of S = 0.29(1), i] = 0.00(1), 
and z = 1.150(15). These values indicate that the bicritical behavior falls in the BAWe 
universality class, for which 6 = 0.285(2), r] = 0.000(1), and z = 1.141(2) 0. 

For Pc < Pc along the A-B coexistence line, a similar analysis shows a crossover from 
the bicritical behavior to sub-critical behavior corresponding to the well known problem 
of the T = one-dimensional kinetic Ising model for which dynamic exponents 6 = 0.5, 
rj = —0.5, and z = 1 are known exactly [^. The two species version of our model, which 
occurs for pc = on the edge of the phase diagram, can be mapped onto this kinetic Ising 
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Model. However, as can be seen in Fig. |, for < pc < Pc a.t short times the dynamic 
critical behavior tends to act more like the bicritical behavior before changing to kinetic 
Ising model behavior at long times. The time which this crossover occurs increases as pc 
approaches p^, but for all pc < Pc the long time dynamical critical behavior corresponds to 
the kinetic Ising model. 



C. Bicritical interface dynamics 

To further analyze the importance of competition in the growth of two equivalent sat- 
urated phases at the bicritical point we also studied the dynamics of an interface between 
those two phases. Starting with a single vacancy between the two domains, we used two 
different methods to analyze the behavior of the interface. Since there must always be at 
least one vacancy between two different saturated phases, in the first method we ignore the 
survival probability P{t) and take 6 = 0. We then measure the number of vacancies in the 
interface {n{t)) oc f and average size of the interface {R'^{t)) oc t^. From 5 x 10^ independent 
runs at the bicritical point, each lasting 10^ time steps, we found the other exponents to be 
rj = 0.285(10) and z = 1.14(2). This type of interface dynamics has been used to study the 
properties of critical interfaces in other models in the BAWe class, where similar results for 



r] and z were obtained ||T8|,|19 



In the second type of interface dynamics simulations, which has not been studied before, 
the simulation is stopped if the interface between the domains has "collapsed" back to one 
vacant site. We introduce a probability of avoiding a collapse P(t) oc t^^ and corresponding 
vacancy concentrations {n(t)) oc f and (i?^(t)) oc . Figure ^ shows results from 10^ 
independent runs each lasting up to 10^ time steps. We find values of 6' = 0.73(2), rj' = 
-0.43(2) and z' = 1.15(2). 

Note the value of the dynamic exponent z or z', which measures the size of the active 
region during surviving runs, is the same in both types of interface dynamics simulations as 
that measured for the defect dynamics. Furthermore, although the exponents S and r] are 
different in the three cases, their sum 6 + 1] (or 6' + rj'), which governs the time evolution 
of the number of vacancies in just the surviving runs, are the same within statistical error. 
This indicates a universal nature of the critical spreading of the active region for models 
with two symmetric adsorbing states which is independent of whether defect or interface 
dynamics is being considered. A similar result holds for some one-dimensional systems with 
infinitely many adsorbing states ||28|| . 

Assuming this conjecture is true, it should be noted that simulations using the first 
type of interface dynamics, where 6 = 0, yield no information beyond that obtainable from 
simulations employing defect dynamics. However, simulations using the second type of 
interface dynamics measure an independent dynamic exponent 6' which we expect to be a 



universal number. Recent measurements on similar models support this conjecture |29 . 



D. Crossover from bicritical to critical behavior 



Finally, we measured the crossover exponent from bicritical to critical behavior. Near 
the bicritical point where the A and B poisoned phases meet, the boundary of the reactive 
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region is expected to behave as {pab — 0.5) oc {pc — PcY, where is the crossover exponent 
2^ . We used the dynamical simulation method to accurately determine the location of 



the DP phase boundary between the reactive phase and the A saturated phase near the 
bicritical point. From the log-log plot of pab — 0.5 versus pc — Pc shown in Fig. |10|, we 
find = 2.1 ± 0.1. Our determination of (p is not as accurate as the other exponents due 
to complications arising from crossover effects. Similar to the crossover from bicritical to 
sub-critical behavior described above, near the bicritical point at short times the dynamical 
behavior is controlled by the bicritical point before changing to the directed percolation 
critical behavior at long times. The time with which the crossover occurs increases as the 
bicritical point is approached, making studies very close to the bicritical point too time- 
consuming. 



V. SUMMARY 

We have studied a simple three species monomer-monomer reaction model to investigate 
the role of symmetry in adsorbing phase transitions. We have shown that, unlike the two 
species monomer-monomer model or the monomer-dimer ZGB model, this model has a re- 
active steady state in one dimension. There are also poisoned states for which the lattice is 
covered by one of the monomer species. The continuous phase transitions between the reac- 
tive phase and the poisoned states meet at bicritical points. Along the first-order coexistance 
line between two absorbing phases the model reduces to the two species monomer-monomer 
model. 

We also constructed a mean-field theory of the model. An unusual feature of the mean- 
field analysis is that the bicritical points lie on the edge of the phase diagram if the corre- 
lations of triplets of adjacent sites are not exactly treated. Only when correlations up to 
triplets of adjacent sites are included does the bicritical point appear inside the phase dia- 
gram, indicating the importance of reproducing the correlations induced by large domains 
of a single saturated phase. 

The dynamic critical behavior at the transition between the reactive phase and a poisoned 
phase is in the DP universality class. At the bicritical points, where there are two equivalent 
poisoned states, the dynamic critical behavior is in the BAWe class. Thus, the universality 
class of the transition changes from DP to BAWe when the symmetry of the adsorbing state 
is increased from one to two equivalent noiseless states. Furthermore, we have shown that 
having a two-fold symmetry in the adsorbing states introduces additional features in the 
dynamics over a model with a unique adsorbing state. In particular, the critical dynamics of 
the interfaces between two different adsorbing states shows a sensitivity to how the dynamics 
is defined, and the survival probability of fluctuations in the size of the interface from its 
smallest value is described by a new universal exponent 6'. However, the critical spreading of 
the reactive region, be it a defect in a single phase or a domain wall between phases, appears 
to be insensitive to the choice of initial conditions. This appears to result from the fact that 
large reactive regions are insensitive to whether the reactive regions are bounded by the same 
or different saturated phases. We do not expect this result to be true in higher dimensions 
where the entropy of domain walls can play a role and nonuniversal critical spreading has 
been observed in other models [pO|] . 
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APPENDIX: RATE EQUATIONS IN THE TRIPLET APPROXIMATION 

The triplet approximation replaces the actual lattice configuration with the average con- 
figuration of each cluster consisting of a three adjacent sites. Define the average number of 
the different number of triplets as 

where is the total number of triplets, which in one dimension is equal to the the number 
of sites, and Nijk is the number of triplets consisting of i, j, and k type monomers {A, B, or 
C) or vacancies {V). The densities of asymmetric triplets, i.e. ijk type triplets with i k, 
are by symmetry assumed to be equal, and are added together. 

In the adsorption controlled limit, triplets with adjacent dissimilar monomers, e.g. A- 
B-B, A-C-V, . . . , are forbidden, leaving 19 allowed types of triplets. However, the triplet 
densities must satisfy 4 separate constraints, which reduce the number of independent triplet 
densities to 15. The first of these constraints, similar to the constraints on the site and bond 
densities given by Eq. (|l]) and Eq. (|) respectively, merely conserves the total triplet density 

ijk 

The other three constraints have no analogues in the site or pair approximations. Because 
each particular lattice site contributes to three different triplets, and the middle and end 
positions of the triplets are not symmetric, the total density of A type monomers occurring 
in say the left position of the triplets must be equal to the the total density of A type 
monomers occurring in the middle position of the triplets 

xavv + '^xava + xavb + xavc = xaav + 2xyAy, 

and similarly for B and C type monomers. 

The equations of motion of the triplet densities can be written as 



where a refers to the event type, and Ax^-J^ are the triplet density changes with an event of 
type a. The different types of events were enumerated above in the discussion of the pair 
approximation. The triplet density changes due to A | and A [ AB ] events are listed in 
Tables II and III, respectively, where 

1 

Vab = xbvv + xbvb + xavb + -xbvc 

and Pijk/xjk is the conditional probability for an i type monomer or vacancy to occur next 
to a jk pair. For example, 

xvvv 

J^vvv/xvv 



Xvvv + k^AVV + ^^BVV + ^^CVV 
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and 



VAV/XAV — 



XyAV 



XvAV + 2^ AAV 

Then taking xyvV) xavv, xbvv, and xcvv to be the dependent triplet densities, the 
equations of motion of the independent triplet densities are 



dx 



AAA 



dt 



Pa xava + C^xava + xavv) Paav/xav 
— {pb Vba + PC Vca) Paav/xav Paaa/xaa 



dXA_ 



av 



dt 



PA 



Xavv + {'^xava + xavv) {Pvav/xav — Paav/xav^ 



+ (PByBA +PcycA) Paav/xav {— ^ + Paaa/xaa — Paav/xaa) 



dx 



VAV 



dt 



Pa xvvv — {'^xava + xavv) Pvav/xav 

— {pb Vba + Pc Vca) {Pvav/xav — Paav/xav Paav/xaa) 



dXAVA 

dt 



Pa —Xava + {'^xwv + xavv) Pavv/xvv 

- {pb + Pc) Xava 

— {pb Vba + Pc Vca) Pvav/xav Pava/avx 



dXA 



VB 



dt 



Pa {'^xvvv + Xavv) Pbvv/xvv 

+PB {'^Xyyv + Xbvv) PaVV/XVV 

—pA xavb — Pa Has Pvbv/xbv Pavb/bvx 
—PB Xavb — Pb Vba Pvav/xav Pavb/avx 
-Pc Xavb 

—Pc ycA Pvav/xav Pavb/avx — Pc Vcb Pvbv/xbv Pavb/bvx 



and similarly for the remaining densities. 
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TABLES 



event type 
Axyy 

^XAA 
^XBB 
/S.XCC 
^XAV 
AXBV 

Axcv 



-2 PA ZAV xvv 
PA ZAV XAV 




Pa ZAvi^xvv - XAv) 






Aj AB] 

hp A XBV (1 + ZAv) 



1 + £Bii 



ipAXBv{l + ZAv)^ 









-^PA XBvi'i- + ZAv)-, 




BV 



TABLE L Bond density changes for different events in the pair approximation. 



event type 



-PA XVVV + (2xyyy + Xavv)PvVV/XVV 
PA XAV A + {2XAVA + Xavv)PaAV/XAV 
XAVV + {2XAVA + XAVv){PvAV/XAV - PaAV/XAv) 
PA XVVV — {2XAVA + XAVv)PvAV/XAV 

PA —XAV A + {^xvvv + xavv)Pavv/xvv 
-XAVV + C^xvvv + XAvv){Pvvv/xvv - Paw/ xvv) 
-PAi^xvvv + xavv)Pbvv/xvv 

-pAi^XvvV + XAVv)PcVV/XVV 

Pa{2xvvv + xavv)Pbvv/xvv 
Pa{2xvvv + XAvv)Pcvv/xvv 



/S.XVVV 

AXAAA 
/S.XAAV 
/S.XVAV 
AXAVA 
AXAVV 
AXBVV 

AXAVB 
AXAVC 



PA 



PA 



TABLE IL Triplet density changes due to an A monomer adsorbing, and remaining on the 
lattice, in the triplet approximation. 
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event type A | AB | 

^XVVV PA XBVV + PA VAB PvBV/XBvi'^ + PbVV/BVx) 

^xbbb -PA VAB Pbbv/xbv Pbbb/xbb 

^xbbv pa VAB Pbbv/xbv (-1 + Pbbb/xbb - Pbbv/xbb 

^XVBV —PA VAB {PvBV/XBV - PbBV/XBV PbBV/XBb) 

^XbVB —PA XbVB - PA VAB PvBV/XBV PbVB/BVX 

^XAVV PA XAVB + PA VAB PvBV/XBV PavB/BVX 

^XBVV PA [xBVB - XBVv] + PA VAB PbBV/XBV 

+PA VAB Pvbv/xbv {Pbvb/bvx — Pbvv/bvx) 

^XCVV PA \xBVC + PA VAB PvBV/XBV PbVC/BVX 

^XAVB -pA XAVB - PA VAB PvBV/XBV PavB/BVX 

^XbvC -PA \xbvC - PA VAB PvBV/XBV PbVC/BVX 



TABLE III. Nonvanishing triplet density changes due to an A monomer adsorbing, and reacting 
with a B monomer to form a AB molecule, in the triplet approximation. 
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FIGURES 




FIG. 1. Phase diagram showing three saturated phases (indicated by the letters), and a 
reactive phase (the unlabeled center region). Solid lines indicate continuous transitions. Dashed 
lines indicate first-order transitions. Bicritical points (filled circles) occur where two critical lines 
meet a first-order line. 
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FIG. 2. Phase diagram in the site approximation. Transitions between the reactive phase 
(unlabeled) and the three saturated phases (indicated by the letters) are continuous. Note that 
unlike the actual phase diagram shown in Fig. |l| the continuous transition lines are straight and the 
bicritical points where two continuous transition lines meet are on the edge of the phase diagram. 
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FIG. 3. Phase diagram in the pair approximation. Transitions between the reactive phase 
(unlabeled) and the three saturated phases (indicated by the letters) are continuous. Note that 
unlike the actual phase diagram shown in Fig. |^ the bicritical points are still on the edge of the 
phase diagram. 
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Pab 



FIG. 4. Phase diagram in the triple approximation. Transitions between the reactive phase 

and the three saturated phases are continuous, while transitions between saturated phases are 
first-order. Inset shows a closcup of the phase diagram near the bicritical point at the end of the 
first-order line separating the A and B saturated phases. All of the qualitative features of the 
actual phase diagram are reproduced in this approximation. 
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FIG. 5. Log-log plot showing the average number of vacancies, ny as a function of time near 
the transition from the reactive phase to the C saturated phase at pab = 0.5. Prom top to bottom, 
the 3 curves correspond to pc = 0.395, 0.39575, and 0.3965. The middle curve corresponds to the 
critical point. Note that the critical line is straight while the other lines have curvature. 
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FIG. 6. Effective exponents using Eq. @ with 6 = 5 for the defect dynamics near the critical 
point at pAB = 0.5 on the line where the C poisoned phase meets the reactive phase. From top 
to bottom, the 3 curves in each panel correspond to pc = 0.3955, 0.39575, and 0.3960, with the 
middle curve corresponding to the critical point. 
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FIG. 7. Effective exponents using Eq. (|6|) with 6 = 5 for the defect dynamics near the bicritical 
point where the A and B poisoned phases meet the reactive phase as defined in Eq. (^) with 6 = 5. 
From bottom to top, the 3 curves in each panel correspond to pc = 0.121, 0.122, and 0.123, with 
the middle line corresponding to the bicritical point. 
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FIG. 8. Average number of vacancies for pAB = 0.5, and pc < Vc showing crossover from 

bicritical behavior to sub-critical behavior. Prom bottom to top, the curves correspond to pc = 0.0, 
0.04, 0.08, 0.10, 0.11, and 0.122. The top curve {pc = 0.122) corresponds to the bicritical point. 
All other curves tend toward a slope of —0.5 at large t. 
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FIG. 9. Effective exponents, as in Fig. ^, for the second type of interface dynamics near the 
bicritical point where the A and B poisoned phases meet the reactive phase. 
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FIG. 10. Location of the critical line as a function of distance from the bicritical point. The 
data falls on a line with a slope corresponding to the crossover exponent ^ = 2.1 zt 0.1. 
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